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Аннотация 

Введение. Усиление антропогенного воздействия на водные объекты требует комплексных решений для оценки 
экологических рисков. В литературе описаны этапы оценки, возможности экологического менеджмента и 
экспертного анализа. Исследуется моделирование рисков в данной сфере. При этом не изучен потенциал 
прогнозирования рисков для состояния воды и биоразнообразия в ходе часто выполняемых гидротехнических 
работ, например дноуглубления. Очевидны актуальность и важное прикладное значение такого подхода. Цель 
исследования — описание математической модели и программного комплекса, которые позволят оценивать 
риски для видового разнообразия экосистемы мелководного водоема при проведении работ в акватории. 
Материалы и методы. Исходной точкой моделирования было описание процесса движения водных масс по 
уравнениям Навье — Стокса и неразрывности при переменной плотности. Уравнение диффузии-конвекции 
использовали для прогнозирования переноса взвешенных и растворенных частиц и оценки влияния примесей 
при эвтрофикации. Создавая алгоритм, задействовали термины и определения, принятые Росстандартом для 
управления рисками в чрезвычайной ситуации. 

Результаты исследования. Для тестирования решения взяли данные о гидромеханических работах в районе 
порта Архангельска. Визуализировали поля концентрации взвешенных частиц, через 0, 15, 30 и 45 минут после 
выгрузки грунта. Установлено, что в процессе оседания взвеси область ее распространения значительно 
расширяется, и это полностью согласуется с данными натурных экспериментов при проведении 
дноуглубительных работ. Рассчитали и свели в таблицу объемы загрязненной воды при отвалах грунта на трех 
участках (при однократном сбросе и в сумме). Для оценки рисков Азовского моря брали максимальные 
концентрации загрязняющего вещества (меди), полученные в ходе замеров, моделирования и дистанционного 
зондирования Земли. В тестах для определения потенциальной опасности вещества исходили из того, что его 
концентрация вызывает реакцию у 50 % организмов. Для рыб потенциально опасная концентрация — 4мг/л при 
длительности влияния 96 ч. Для зоопланктона — 50 мг/л и 48 ч. Для микроводорослей 20 мг/л и 72 ч. Получено 
значение нормализованного риска — К, = 0,52. Признан значимым риск концентрации меди 80 мкг/л в водах 
Азовского моря. Выявлена тенденция увеличения солености Азовского моря и стратификация водных масс по 
содержанию кислорода, что согласуется с результатами экспедиционных исследований. 

Обсуждение и заключение. Разработанный подход позволил оценить изменение качества вод Азовского моря и 
описать некоторые трансформации акватории. Речь идет, в частности, о распространении взвешенных частиц и 
о районах их оседания. Указанные процессы обусловливают изменение рельефа дна, которое, в свою очередь, 
может привести к сокращению видового состава водоема. 
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Введение. Развитие экономики сопровождается масштабными работами на территории водных объектов. 
В качестве примеров можно привести строительство Золотого, Русского и Крымского мостов, дноуглубление [9 
целью расширения акватории для судоходства, прокладку подводных Газо- и нефтепроводов и многое другое. 
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Обязательная часть таких проектов — оценка потенциального воздействия на экологию, прогнозирование 
рисков. В России эти вопросы регулируются Требованиями к материалам оценки воздействия на окружающую 
среду'. Кроме того, действует ГОСТ Р ИСО 14001-2016. Это официальный перевод международного стандарта 
[ЗО 14001:2015, определяющего деятельность предприятия и его системы экологического менеджмента при 
формировании заявки и реализации проекта. Стоит отметить, что действующая редакция стандарта направлена 
в большей степени на предупреждение чрезвычайной ситуации, а не на ликвидацию ее последствий. 

Авторы работы [1] выделяют три ступени оценки риска для окружающей среды, в [2] описаны пять этапов. 
Исследование [3] фокусируется на рисках нефтяных и газовых компаний. Здесь же рассматриваются 
позитивные аспекты внедрения системы экологического менеджмента. 

Значительная часть материалов посвящена развитию систем корпоративного экологического менеджмента 
для обеспечения международной или национальной безопасности в сфере экологии [4]. В[5] изучена 
корреляция экологической ответственности компании и стоимости ее акций на фондовом рынке. Аналогичная 
проблематика рассматривается в [6]. Речь идет о том, что экологическая ответственность французской 
компании сказалась на росте ее рыночной стоимости. Оценивалось влияние на капитализацию премии «Грин 
компани эвордс» (Стееп Сотрапу А\маг4$), которой отмечают лучшие с точки зрения экологии объекты и 
предприятия [7]. Во многих случаях учитывается экспертный подход к оценке экологических рисков 
компаний [8]. Однако, у него есть существенные недостатки, связанные с субъективностью выводов и 
сложностью обработки большого объема информации в условиях неопределенности [9]. Поэтому более 
предпочтительны математические методы оценки риска. 

В [19] описываются популярные подходы к оценке рисков с акцентом на индексный. В работе [11] рассмотрена 
модель для оценки экологических рисков с использованием теории доказательств Демпстера — Шафера, которая 
хорошо себя зарекомендовала при решении задач в условиях неопределенности. Тем не менее ее практическое 
применение осложняется необходимостью многих сложных расчетов, в том числе комбинированных. В статье [12] 
описана математическая модель на основе вероятностного подхода и определения интегральной величины оценки 
риска. В литературе не рассматриваются возможности прогнозирования комплекса рисков для состояния воды и 
биоразнообразия при выполнении гидротехнических работ. Однако именно такой подход следует признать 
актуальным, имеющим очевидное прикладное значение на современном этапе развития экономики. 

Цель представленного исследования — описание математической модели и программного комплекса, 
которые позволяют оценивать уязвимости и риски для гидрохимических параметров акватории и видового 
разнообразия экосистемы водоема при углублении дна мелководного водоема. В первую очередь речь идет о 
прогнозировании распространения и осаждения взвешенных частиц. Эти процессы существенно влияют на 
рельеф донной поверхности. Кроме того, моделируются зоны распространения загрязняющих веществ. Их 
токсичность может стать фактором сокращения видового разнообразия водной экосистемы. 

Материалы и методы. Для прогнозирования результатов антропогенного влияния на экосистему водоема 
(например, при проведении дноуглубительных работ) предлагается комплексная математическая модель 
движения взвешенных частиц в водной среде. Она учитывает ветровые течения, движение водной среды и 
стоки рек, геометрию акватории, турбулентный обмен, переменную плотность водной среды, которая может 
зависеть от солености или взвеси. При моделировании предлагается для каждой фракции взвеси принимать 
во внимание скорость осаждения, определяемую размером и формой частиц. Решение можно дополнить 
моделью эвтрофикации водоема, позволяющей учитывать влияние биогенных веществ на процессы, 
протекающие в водоеме. 

Постановка задачи 

Модель движения вод. Для описания процесса движения водных масс в акватории задействуем модель 
гидродинамики [13]. Она включает в себя представленные ниже выражения. 

1. Уравнения движения Навье — Стокса: 

ди ди би ди _ ТОР °( Э > м "(| м. 
и +у У + | + | + У : 
01 Ох ду 02 р ах дх\ дх ду д 02\ 02 
ду ду д ду ТОР Од[ д\ 09| 9] д[| д» 
и чу и + | )* |» } (. } (1) 
01 Ох ду 02 роду дх\ дх] д»\ д) 02\ 02 
ди ди 0» 0» _ ТОР 2 ( =) "(| ==] (| =) 
и У + + | + |п | У + 2. 
01 Ох ду 02 р 02 дх\ дх ду д» 02\ 02 


1 Об утверждении требований к материалам оценки воздействия на окружающую среду. Приказ Министерства природных ресурсов и 
экологии Российской Федерации. ЧВГ: Врз://40с$.сп.ги/Чоситаен/ 573339130 (дата обращения: 26.02.2024). 
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Здесь /= {и, у, \} — вектор скорости движения водной среды [м/с]; Р — давление [Па]; р — плотность [кг/м]; 
и, у — горизонтальная и вертикальная составляющие коэффициента турбулентного обмена [м?/<]; 
$ — ускорение силы тяжести [м/с?]. 

2. Уравнение неразрывности в случае переменной плотности: 


др, (ри) д(ру) О(ру) _ 
01 дх ду 02 


Начальные и граничные условия для системы уравнений (1) и (2) подробно описаны в [13]. 
Модель распространения взвешенных частиц. Для прогнозирования переноса взвешенных и 
растворенных частиц используем уравнение диффузии-конвекции [14]: 


дс, д(ис,) д(у,) 9((»+ у.» )е,) 0[ дс, 0[ бе, 9[ 0с, 
+ + + = +—|ы +—у 
Ге дх ду 02 0х\ Обх] д. 09] 42\ 02 
где с, — концентрация 7-й фракции взвеси [мг/л]; и’, — скорость гравитационного осаждения 7-й фракции 
взвеси [м/с]; РЁ’, — функция интенсивности источников 7-й фракции взвеси [мг/(л-с)]. 


Система уравнений (3) рассматривается при начальных и граничных условиях, подробно описанных в [14]. 
На основе системы (1)-(3) можно смоделировать процессы движения и осаждения взвешенных частиц при 


0. (2) 


Е,, (3) 


проведении дноуглубительных работ, а также рассмотреть возможность оптимизации площадей отвалов 
грунта. Модели целесообразно задействовать для сокращения вреда экосистеме водоема. 

Модель эвтрофикации. Для оценки влияния примесей рассмотрим модель эвтрофикации вод на основе 
уравнения диффузии-конвекции [15]: 


д 5 
85, ‚ 0(иб.) 9(5,) д((*+м»,)5.) 0[ 88,\, 9[ 25, 9[85, 


+Р,., (4) 


и м д р т а и 


где 5, — концентрация 7-й примеси [мг/л]; у,, — скорость гравитационного осаждения 7-й примеси [м/с]; 
Е, — функция, описывающая интенсивность химико-биологического источника 7-й примеси [мг/(л`с)]. 

Система уравнений (4) рассматривается при начальных и граничных условиях, которые приводятся в [15]. 

Можно использовать (4), в частности, для описания: 

— влияния примесей, поступающих в водоем в результате антропогенного воздействия, например при 
аварийной ситуации или сбросе сточных вод; 

— процессов окисления и восстановления марганца; 

— окисления сероводорода; 

— расчета концентрации растворенного в воде кислорода. 

Данные о размещении источников можно получить с помощью дистанционного зондирования Земли. 

Алгоритм оценки риска при антропогенном воздействии. Известные подходы к пониманию риска 
систематизированы в [16]. В данной работе используется определение риска в соответствии с 
ГОСТР 55059-20122. Риск — это вероятность наступления чрезвычайной ситуации или ожидаемый размер 
ущерба, связанный с реализацией неблагоприятного исхода. Рассмотрим пример, когда в воду поступают 
примеси, негативно влияющие на экосистему: 

В=Е/у Е. (5) 

Здесь К — фактор риска; Ех — потенциальная опасность или токсичность примеси; Ез — фактор 
экспозиции или количество примеси в перерасчете на одну биологическую мишень. Для водоема фактор 
экспозиции можно рассчитать по формуле: 

Е; =5-Ру-Ёвл, (6) 
где 5 — заданная или рассчитанная на основе модели (4) концентрация примеси; Ру — индекс скорости 
гидролиза; Авх — фактор биоаккумулирования? или накопления (например, тяжелых металлов) [17]. 

Для определения потенциальной опасности примеси или загрязняющего вещества Ек в мировой практике 
принято использовать набор стандартизированных тестов“. В странах Европейского Союза выделяют 
описанные ниже количественные оценки токсичности вещества. 


2 ГОСТР 55059-2012. Безопасность в чрезвычайных ситуациях. Менеджмент риска чрезвычайной ситуации. Термины и определения. 
Москва: Стандартинформ, 2018. 8 с. 
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— МОЕС — по оБзегуе4 еЁесф сопсепгайоп, максимальная недействующая концентрация вещества. Как 
правило, эта концентрация по сравнению с контрольной не вызывает статистически значимого негативного 
эффекта (вероятность его наступления не превышает 0,05) в течение заданного времени воздействия. 

—ГОЕС — ю\е$ обзегуе4 еЁесё сопсепбайоп, наименьшая наблюдаемая эффективная концентрация. Ее 
действие вызывает ответную реакцию тестируемых организмов. Это могут быть нарушения обмена веществ, 
роста, развития, размножения и даже гибель. Значения ГОЕС выше значений МОЕС. Если известен эффект 
(процент эффекта) ГОЕС, можно примерно определить МОЕС: МОЕС =ГОЕС/2 (при 10 % <ГОЕС < 20%). То 
есть если ГОЕС оказывает влияние на 10—20 % исследуемой популяции, то значение МОЕС можно примерно 
определить как половину от ГОЕС. 

— МАТС — тахипа! ассераЫе ‘юх1сап{ сопсептайоп, максимально допустимая концентрация токсичного 
вещества. Это расчетный параметр, определяемый как среднее геометрическое МОЕС и ГОЕС; 

— ЕС, — еНес{ сопсешгайоп (можно перевести как «концентрация эффекта, при которой наблюдается х % 
эффекта по сравнению с контрольной группой»). То есть это концентрация вещества, при которой наблюдается 
ответная реакция у х% тестируемых организмов. Например, при концентрации ЕС5о реагируют 50% 
организмов. Для расчета данного критерия часто используется статистический метод (например, 
регрессионный анализ). При этом необходимо использовать достаточное количество групп концентраций (доз), 
поскольку точность оценки зависит от числа и интервала концентраций, а не от размера выборки для каждой 
концентрации. Наряду с аббревиатурой ЕС, используют СС, или Г(Е)С,. 

В рамках данного исследования рассматривались три набора тестов“, применимые к природным водным объектам. 

Результаты исследования. На основании рассмотренных моделей и подходов к оценке риска 
сформировали описанный ниже алгоритм из пяти шагов. 

Шаг 1. Изучение характеристик акватории (географических, климатических, гидрологических). На данном 
шаге также рассматриваются потенциальные источники поступления взвеси и примесей. Это могут быть 
работы по расширению и очистке акватории, стоки рек, промышленные сбросы, судоходство ит.д. Следует 
определить возможный объем взвеси. 

Шаг 2. Моделирование сценариев распространения взвесей и примесей на основе уравнений (1)-(4). 

Шаг 3. Оценка токсического эффекта Еф по трем выбранным тестам. 

Шаг 4. Оценка фактора риска для каждого вещества на основе (5)—(6) и расчет нормализованного фактора 
риска Кл Е [0, 1]. Для нормализации используется формула (5 — Зив)/($ах — Эта), ГДе Эта И Эщшах — минимальная 
и максимальная концентрации рассматриваемого вещества. 

Шаг 5. Формулирование предварительных выводов. Будем считать риск высоким при К, > 0,55, значимым 
при 0,3 < К, < 0,55 и отсутствующим при А, < 0,3. 

Шаги 2-4 могут повторяться, если нужно определить оптимальный вариант проведения работ. Например, 
при дноуглублении можно моделировать распространение взвешенных веществ для различных точек сброса 
грунта и оценивать вред, причиняемый экосистеме. 

Для численной реализации моделей (1)-(4) и описанного выше алгоритма оценки риска при антропогенном 
воздействии разработан программный комплекс на языке С++. Он объединяет четыре программных модуля. 
Ниже описано назначение каждого из них. 

1. Расчет трехмерного вектора скорости движения водной среды на основе системы уравнений (1)-—(2). 
Принимается во внимание сложная геометрия расчетной области. 

2. Расчет на основе системы уравнений (3) переноса взвешенных частиц и их оседания. В данном случае 
имеется в виду многокомпонентная взвесь. 

3. Моделирование эвтрофикации водоема на основе системы уравнений (4). Учитываются фито-, 
зоопланктон, а также 13 химических элементов и их соединений, которые существенно влияют на 
гидробиологические процессы водоема. Это растворенный кислород, железо, сера, сероводород, сульфаты, 
сульфиты, азот, аммонийный азот, нитриты, нитраты, фосфаты, силикаты, кремнекислота. 

4. Оценка рисков антропогенного воздействия на основе описанного выше алгоритма. 

Моделирование процессов движения и оседания взвеси при дампинге грунта. Разработанный 
программный комплекс моделирует процессы заиления судоходных каналов Дона. К тому же он позволяет 
изучать транспорт взвешенных частиц и переформатирование донной поверхности Азовского моря. Особое 
внимание уделяется прибрежной зоне и устьевым районам. Кроме того, софт использовался для оценки 
экологического состояния акваторий при проведении дноуглубительных работ в Белом море, в районе 
Двинской губы. 
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В качестве примера рассмотрим моделирование процессов движения и осаждения взвеси при проведении 
работ по расширению акватории [18]. С этой целью задействуем данные по работам в районе порта 
Архангельска. Для моделирования распространения и осаждения взвешенных частиц в качестве исследуемой 
области водоема взяли участок длиной 3 км вдоль направления течения. Его ширина — 1,4 км; глубина — 10 м. 

Физические параметры водной среды и взвеси: 

— скорость течения — 0,2 м/с; 

— плотность взвешенного материала — 1600 кг/мз; 

— скорость осаждения взвеси — 2,042 мм/с; 

— содержание частиц диаметром менее 0,05 мм в грунте — 26,83%; 

— объем выгружаемого насыпного материала — 741 м3. 

Параметры расчетной области: 

— шаг по горизонтальным пространственным координатам — 20 м; 

— шаг по высоте — 1 м; 

— расчетный интервал — 2 часа; 

— шаг по времени — 1 минута. 

На рис. 1 представлены поля концентрации взвешенных частиц (в мг/л), соответствующие разным 
интервалам времени. Трехмерный срез показывает сечение расчетной области плоскостью, которая проходит 
через точку выгрузки и образуется векторами, направленными вертикально и вдоль течения (течение — 
слева на право). 


2, М 
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Рис. 1. Поля концентрации взвешенных частиц (в срезе «глубина — ширина» водоема) в разные интервалы времени с 
момента выгрузки грунта: а — Т= 0; 6 — Т= 15 мин; в — Т= 30 мин; г — Т =45 мин. Цветная легенда показывает 
концентрацию взвеси. На горизонтальной оси — данные о ширине водоема, на вертикальной — о глубине 
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На основе полученных результатов рассчитали объем загрязненной воды при отвалах грунта на трех 
участках сброса (таблица 1). 


Таблица 1 
Объемы загрязненной воды при сбросе грунта на трех участках, млн м3 
Общий объем загрязненной воды 
Объем загрязненной воды при одном 
сбросе с учетом концентрации с учетом концентрации взвеси 
взвеси в воде Обкем 
Участок в воде 
сбросов Всего 
от 0,25 от 20 вол88100 от 0,25 от 20 более 100 
до 20 мг/л |до 100 мг/л мг/л до 20 мг/л | до 100 мг/л мг/л 
1 0,89 0,245 0,15 124 110,36 30,38 18,6 159,34 
2 0,813 0,202 0,105 50 40,65 10,1 5,25 56 
3 0,889 0,24 0,15 45 40,01 10,8 6,75 57,56 


Из таблицы видно, что при сбросе грунта высокая концентрация взвеси фиксируется в относительно малом 
объеме воды (это, допустим, часть воды, которая перемешана с грунтом). Большее количество воды имеет 
более низкую концентрацию (как, например, вода в непосредственной близости от сброшенного грунта). 
Значительные объемы воды загрязняются с меньшей концентрацией взвеси в процессе сброса грунта, его 
осаждения и переноса течением. 

Участки сброса грунта и количество сбросов определялись по техническому заданию для работ по 
углублению акватории порта. Выявлены области водоема, в которых возможна гибель фито- и зоопланктона. 
Рассчитаны площади этих областей. 

Оценка экологического состояния акватории Азовского моря. Экологическое состояние мелководного 
водоема оценивалось на примере Азовского моря. При исследовании использовали два подхода. Первый 
предполагает анализ баз постоянно пополняющихся данных. Это результаты экспедиционных изысканий 
авторов данной статьи в Азово-Черноморском бассейне. Второй связан с применением программного 
комплекса, который моделирует гидродинамику и гидробиологию водоема, а также процессы транспорта 
взвешенных частиц. По результатам сравнения полученных оценок проводились калибровка и верификация 
математических моделей для более точных оценок рисков. 

Программный комплекс дает возможность оценить риски таких явлений, как: 

— превышение предельно допустимых концентраций опасных загрязняющих веществ; 

— бурный рост в процессе эвтрофикации вредоносных и токсичных водорослей («цветение» вод). 

Таким образом, задействовав новый софт, можно сократить затраты на экспедиции, которые организуются 
для отбора проб воды. 

Отметим также, что для контроля экологического состояния изучаемого водного объекта использовались 
данные дистанционного зондирования Земли. 

Приведем пример обработки данных замеров, полученных в экспедиции. Рассмотрим процесс поступления 
меди от предприятий металлургической промышленности со стоком реки Дон в Азовское море. Медь — это 
тяжелый металл. Она может накапливаться в организме до опасных показателей. 

При оценке рисков мелководного водоема (Азовского моря) использовались максимальные значения 
концентрации загрязняющего вещества, полученные в ходе замеров, математического моделирования и данных 
дистанционного зондирования Земли. В фиксированной точке водоема или исследуемой области выбирали 
конкретное измеренное или рассчитанное значение концентрации загрязняющего вещества. Предположим, что 
замеры выявили фактическую концентрацию взвеси в воде 80 мкг/л. В зависимости от участка и глубины 
забора проб концентрация меди в водах Азовского моря колеблется от 0,001 мкг/л до 100 мкг/л. 

Индекс скорости гидролиза для меди примем равным Рх = 1,5, а фактор биоаккумулирования (накопления 
вещества) — Ави = 2. Тогда согласно (6) фактор экспозиции Ё5 = 240 мкг/л. 

Как отмечалось выше, для определения потенциальной опасности вещества Ё/5 рассматриваются три теста: 

— [Со в течение 96 часов для рыб (судак) — 4 мг/л; 

— [Сэ в течение 48 часов для зоопланктона (дафнии) — 50 мг/л; 

— [Сэ в течение 72 часов для ингибирования роста микроводорослей — 20 мг/л. 
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Примем во внимание (5) и описанный выше алгоритм. Получим значение нормализованного риска К, = 0,52. 
Таким образом, риск наличия меди в концентрации 80 мкг/л в водах Азовского моря можно охарактеризовать 
как потенциально значимый. 

Также стоит отметить, что модель (1)-(4) позволила отследить тенденцию увеличения солености Азовского 
моря и стратификацию водных масс по содержанию кислорода. Это согласуется с результатами проведенных 
ранее экспедиционных исследований [19]. 

Остановимся отдельно на реализации 4-го шага алгоритма, то есть на оценке фактора риска для каждого 
загрязняющего вещества. Используя систему (5)—(6), расчет нормализованного фактора риска А, Е [0,1] и 
разработанный программный комплекс, вычислили максимальные значения концентраций основных 
загрязняющих веществ, характерных для Азовского моря. 

Завершающий анализ рисков проводили по итогам обработки экспедиционных данных и результатов 
математического моделирования. Судя по индексу загрязненности воды, экологическое состояние Азовского 
моря улучшается. Ранее вода определялась как «значительно загрязненная», теперь как «умеренно 
загрязненная» [19]. 

Обсуждение и заключение. Предложенное софтверное решение дает возможность прогнозировать 
результаты антропогенного влияния на динамику качества вод Азовского моря. Разработанный программный 
комплекс позволяет моделировать распространение в акватории взвешенных частиц, определять основные 
районы их оседания. Такие процессы меняют рельеф дна. Следствием может быть существенное сокращение 
видового состава в районе оседания взвеси. Согласно полученным результатам, уменьшение площади отвала 
грунта при проведении дноуглубительных работ способно ограничить вред, причиняемый экосистеме водоема. 
Ориентируясь на полученные таким образом данные, можно в целом судить о потенциальных негативных 
последствиях, связанных с экономическим ущербом и рисками для здоровья людей. 
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